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Abstract 

We investigate the problem of the superuniversality of the phase transition 
between different quantum Hall plateaus. We construct a set of models which 
give a qualitative description of this transition in a pure system of interacting 
charged particles. One of the models is manifestly invariant under both Dual- 
ity and Periodic shifts of the statistical angle and, hence, it has a full Modular 
Invariance. We derive the transformation laws for the correlation functions 
under the modular group and use them to derive symmetry constraints for 
the conductances. These allow us to calculate exactly the conductivities at 
the modular fixed points. We show that, at least at the modular fixed points, 
the system is critical. Away from the fixed points, the behavior of the model 
is determined by extra symmetries such as Time Reversal. We speculate that 
if the natural connection between spin and statistics holds, the model may 
exhibit an effective analyticity at low energies. In this case, the conductance 
is completely determined by its behavior under modular transformations. 
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I. INTRODUCTION 



The problem of the nature of the phase transitions between the plateaus of the Quantum 
Hall States is one of the most challenging questions in the physics of strongly correleted Fermi 
systems. These phase transitions separate states whose most important physical feature is 
their inherent topological nature. Thus, it has been shown that the fluid ground states of 
Quantum Hall systems have a form of topological order, instead of conventional long range 
order |JJ. Furthermore [Q-^J], the low- lying excitations are fractionally charged quasiparticles 
with fractional statistics 0. A direct consequence of the topological order is the universal 
properties of the edge states f|. In addition to the values of the Hall conductance, the 
plateaus are uniquely defined by the quantum numbers of their excitations. Hence, the 
effective theories which describe the plateau behaviors must be topological field theories. It 
has been known for some time that Chern-Simons gauge theories give an accurate description 
of the plateaus . 

It is then natural to ask if the behavior of the system at the phase transtition between a 
pair of plateaus depends on the properties of the plateaus. A theory of the transition between 
Hall plateaus was proposed recently by Kivelson, Lee and Zhang (KLZ) ||. This theory 
explains the rich phase diagram seen experimentally. KLZ used a bosonic description of the 
FQHE and, in this picture, the transition between plateaus becomes strongly reminiscent to 
the superfluid-insulator phase transition in a disordered system. The transition beteween 
plateaus appears to be generally a first order transition in the pure system but it is smoothed 
to a second order transtion by effects of disorder. In this picture, the stable FQHE plateaus 

2 

are the phases in which the Hall conductivity a xy = v\ has precisely defined values. Here 
v is a fraction whose denominator is always odd. In this state the longitudinal conductivity 
vanishes exactly, a xx = . Since in such theories the Hall conductivity is the coefficient of the 
effective Chern-Simons term for the external electromagnetic field, the transition between 
plateaus becomes a transition between states with effective actions for the electromagnetic 
field with different Chern-Simons coefficients. At the transition a xx ^ (and perhaps 
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universal) and the denominator of a xy is even. In the language of the renormalization 
group, the plateaus are represented by stable fixed points described by effective actions 
which have a pure Chern-Simons form PJlQ[|. From this viewpoint it is natural to ask if 
the critical properties at the transtion between plateaus depends on the precise values of 
the Chern-Simons coefficients of the fixed points or if there is some sort of superuniversality 
at work. Kivelson, Lee and Zhang || have argued that the localization length exponent 
at the transition between plateaus in a disordered quantum Hall system should be the 
same for all these phase transitions, while the longitudinal conductivity at the transition 
may be universal although different for different transitions. Hence, KLZ have argued that 



these transtions are superuniversal. Furthermore, Liitken and Ross [fTT]JT^1 have conjectured 
that the theories for the FQHE posses an effective modular invariance which seemingly 
implies such a superuniversality in the sense that all the critical points have the same 
critical exponents and the critical conductances are related by modular transformations. 
Recent Monte Carlo simulations on a related system by Manousakis ]13[ seem to support 
the notion of superuniversality. 



Recent experimental work by Shahar et. al. [H] on the transition between quantum Hall 
and insulating states have revealed an effective self-duality relating the properties of the 
system on either side of the phase transition. This self-duality is observed in a suprisingly 
broad neighborhood of the critical point. This observation gives strong motivation to study 
systems with particle- vortex duality. Of course such a symmetry can emerge near the phase 
transition even if the microscopic theory is not self-dual. 

Chern-Simons gauge field theories thus provide a natural description of the plateaus of 
the fractional quantum Hall effect (FQHE) as well as of other problems in condensed matter 
physics and field theories in two space dimensions. Chern-Simons theories appear in the 



effective theories for the low energy physics R|T5|,|T6[ (in the Landau-Ginzburg sense). Their 



presence makes manifest the breaking of time reversal invariance in the underlying system. 
Chern-Simons theories are also the natural framework for statistical transmutation JT7|-pT 



In particular, in the FQHE Chern-Simon theories appear in the form of exact mappings 



between theories in which the particles which bear the same electric charge but, in general, 
have different statistics ]7|,|22],[23| . This latter approach has led to two mutually complemen- 
tary descriptions of the FQHE and its hierarchies. One is based on a bosonic picture within 
which a superfluid analogy plays a central role. The other uses a fermionic picture with an 
integer Hall effect as the reference analogy |M . 



Motivated by these considerations, in this paper we consider the critical properties of 
phase transitions in theories of matter coupled to a Chern-Simons gauge field. Our main 
goal is to investigate if the critical properties of such systems are altered in any relevant 
way by the coupling of the order parameter field to a fluctuating Chern-Simons gauge field. 
In contrast with the generally accepted assumption that the transition between plateaus 
is necessarily a first order transition in a pure system, we construct a model in which the 
transition is apparently second order, even in the abscence of disorder. This model can serve 
as the prototype of the phase transition between plateaus. The strategy that we follow in this 
work is to construct physically reasonable models which manifestly have these symmetries 
and to use them to describe the plateaus and the phase transitions. The models should also 
yield the correct electromagnetic response functions at its stable fixed points which will then 
be identified with the plateaus. 

More specifically, we consider systems of interacting particles with fractional statistics 
parametrized by an angle 5. In Chern-Simons language, this corresponds to a coupling 
constant 9 = 1/25. In order to address the issue of superuniversality it is necessary to define 
carefully the symmetries exhibited by these systems, even at the microscopic level. The 
problem discussed here is the same as to ask if the universality classes of a phase transition 



depends on the fractional statistics f25fl . The notion of fractional statistics is clearly well 
defined for systems with a finite number of particles (first quantization). In this case, it 
is straightforward to check that all theories of this type are invariant under a shift of the 
statistical angle 5^5 + 2nk (where k is an integer) or, equivalently, -s — > g + 4:iik. This 
is the invariance of the theory of the FQHE under the attachment of an even number 2k 
of flux quanta to each particle. We will call this symmetry Periodicity and denote it by P 



(not to be confused with parity!). Also, in systems in which particle-hole symmetry is exact, 
the theory is invariant under a change in the sign of the Chern-Simons coupling constant, 
9 — ► —9. For obvious reasons, we will call this symmetry Charge Conjugation and denote 
it by C. Finally, we will show that there is another symmetry at work in these systems, 
Duality, which we will call S. However, as in other situations in statistical physics, Duality 
is not a symmetry but rather it maps different theories among each other. We show however 
that it is possible to construct exactly self-dual theories which will serve as fixed points. 

In this paper we construct two classes of models, which we call Model I and II respectively. 
Both models exhibit quantum critical behavior at zero temperature, In Model I, at least, 
this behavior is associated with zero temperature phase transitions. Model II appears to be 
critical at zero temperature for all physically meaningful values of the coupling constants. 
Both types of model are explicitly Periodic in the statistical angle and, in addition, Model 
II is explicitly self-dual. 

Model I describes a system of charged particles with (in general) fractional statistics 
and short range interactions. These models are not self-dual in the sense that the form of 
the interaction is not in general preserved by the Duality transformation. In fact, there 
are cases in which a theory of particles with short range interactions is dual to a theory of 
vortices (the particles of the dual theory) which have long range interactions. This is in fact 
the familiar case of a superfluid-insulator transition. In general, the stable phases of the 
systems described by Model I have the qualitative properies of Quantum Hall plateaus and 
it is qualitatively similar to the Landau- Ginzburg theory of the FQHE. This model gives the 
correct plateau behavior. However, in general, it is not self-dual and its properties at the 



phase transtion betweeen plateaus are not understood. Recently, Pryadko and Zhang |26 
and Balachandran, Chandar and Sathialpan |27| have considered a model closely related to 
Model I and discussed its properties under duality transformations. Our conclusions on the 
physics of these model agrees with theirs. 

The physics of the transitions between Quantum Hall plateaus strongly suggests that, at 
least very close to the transition, the system should be effectively self-dual. This motivates 



the construction of a model, Model II, which is explicitly self-dual. However, Model II will 
turn out to have interactions which are long ranged and retarded. Physically, Model II 
describes the behavior of a system whose excitations are charged particles with fractional 
statistics, restricted to move on a plane. The long-range retarded interactions are due to 
the quantum fluactuations of the electromagnetic field in the full 3 + 1-dimensional space- 
time [ 28] . Naturally, when the coupling constant of the long-range interactions becomes very 



small, the short range interactions (of the sort described by Model I) become dominant. 

We will show that Model II exhibits quantum critical behavior of a sort that is readily 
analyzed. We conjecture that this behavior is characteristic of the behavior of a broader 
class of systems at their zero temperature critical points, including Model I. Model II is 
manifestly periodic under shifts of the statistical angle 5. We find it convenient to define 
a parameter / = 5/2ir. The interactions between these particles is represented by a long- 
ranged kernel which at long wavelenths we choose to be of the form 2n g / ykp, where is 
the three-momentum (i. e., momentum and frequency). The parameter g is be the coupling 
constant and, for stability, it must be positive (repulsive) . This model is both self-dual and 
periodic. We show that the partition function of Model II is invariant under the Duality 
transformations 

9 -> 2 ; f2 /-> — 3-7-72" (1-1) 

g 2 + P t + P 

as well as under Periodicity 

/->/ + ! (1-2) 

These two transformations do not commute with each other. Instead, they generate an 
infinite discrete group known as the Modular Group or SL(2,Z). Thus these systems are 
invariant under that action of SL(2, Z). In fact, in terms of the complex variable z = f + ig 
(which has Imz = g > 0), the group of Modular transformations becomes the group of 
fractional linear transformations of the complex plane with integer coefficients 

az + b 



cz + d 
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;i.3) 



where a, b, c and d are integers such that ad — be = 1. 

Since Model II is manifestly invariant under both Periodicity and Duality, its partition 
function is Modular Invariant. The Modular Group SL(2,Z) has a rich structure. In 
particular it has an infinite set of fixed points, namely points of the upper half plane which 
are invariant under at least one, non-trivial modular transformation. We will refer to these 
points as the modular fixed points. Liitken and Ross fl2| have used this structure as the 
basis for a conjecture concerning the behavior of possible renormalization group flows in 
QHE systems. However, in the work of Liitken and Ross, SL(2, Z) acts directly on the 



complex conductance a xy + io~ xx . Consequently, the fixed points of SL(2, Z) in reference (12 
are the allowed values of the complex conductances. In contrast, in Model II the modular 
fixed points are points of special symmetry. In analogy with the Kramers- Wannier argument 



for the two-dimensional Ising Model |29j , we can expect the phase transitions to occur at the 



modular fixed points. Nevertheless, we will show below that Modular Invariance is enough 
to determine the conductance of Model II, at least at the modular fixed points. However, 
we will also find that there are other values of / and g where the system is critical. 

Having an explicit construction of a model which is modular invariant has allowed us to 
investigate the consequences of modular invariance on physical observables. In particular, 
we show that the current correlation functions of these models transform (almost) like a 
modular form. More precisely, we show that from the current- current correlation function of 
this system it is possible to construct an amplitude, which we call D(z), which determines 
the complex conductance. By construction, D(z) is Periodic under z — > z + 1 and, under 
Duality, it transforms like D(—l/z) = z 2 D(z) + z. Combining both transformations we find 
that, under a general modular transformation, D(z) transforms like 

D(^±^) = (cz + d) 2 D(z) + c(cz + d) (1.4) 

This functional relation is a direct consequence of the modular invariance of the partition 
function and of the fact that this invarience is broken by an external electromagnetic field in 
a well defined manner. In fact, although we derived this functional relation in the context 



of Model II, we expect it to be generally valid for any modular invariant system. What is 
very interesting is the fact that this functional relation is sufficient to determine the val- 
ues of the conductances for all the fixed points of the modular group SL(2, Z). Below, we 
derive a formula which yields the allowed values of the conductance at the modular fixed 
points. Furthermore, by making additional assumtions on the possible analytic properties of 
the correlation functions compatible with the symmetries of the system, we can also deter- 
mine their behavior away from the modular fixed points. In particular, modular invariance 
requires that, quite generally, a xx is finite, at least at the modular fixed points. 

Interacting charged particles without impurities are generally expected to be either insu- 
lators or perfect metals (or even superconductors). Correspondingly a xx = for an insulator 
while a xx = oo for a perfect metal (or superconductor). In order to find a xx 7^ (at zero 
temperature) it is usually assumed that the system must be disordered (or break galilean 
invariance in some way). However, there is an alternative scenario which yields a finite 
a xx in a non-disordered system: a Quantum Critical Point. At a Quantum Critical Point, 
the spectral functions of all the correlation functions are changed drastically by the effects 
of fluctuations on all length scales. In particular, the pole structure associated with the 
existence of well defined elementary excitations is wiped out and replaced by branch cuts. 
Hence, the spectrum becomes a broad incoherent continuum. A direct consequence of this 
behavior is that a xx can now be finite. Thus, we will interpret having a finite o xx in a (pure) 
system of charged interacting particled as a signature of a Quantum Critical Point. 

Some time ago C. J. Callan and D. E. Freed found in their work on the Dissipative 



Hofstadter Model |30| an algebraic structure very similar to the one we use in this paper. It 
is clear that Model II is quite analogous to the Dissipative Hofstadter Model. In particular, 
our coupling constant g plays the role of the dissipation in the Dissipative Hofstadter Model. 
Similarly, the Modular Group also played a central role in the work of Shapere and Wilczek 
]J and of Rey and Zee [[HJ on self-dual Zn models in 2 + 1 dimensions. 
The paper is organized as follows. In Section |TJ we introduce a class of models which 

we 



are manifestly invariant under periodic shifts of the statistical angle. In Section [III 
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construct the Duality transformation. In particular, we discuss the way the presence of 
external electromagnetic perurbations alter the way these systems transform under Duality 
In Section [TV] we discuss the electromagnetic response functions and their properties under 
Duality. In Section [V] we introduce Model I and Model II and discuss the nature of their 
phase diagrams. This section contains the most important results of this paper, including 



the determination of the complex conductance at the modular fixed points. In Section [VI 



we 



discuss the role of periodicity in quantum field theories of matter coupled to Chern-Simons. 



Section VII is devoted to the conclusions. 



II. MODELS 

For all the reasons invoked above we will adopt a lattice formulation in which all symme- 
tries are manifest. We will work in discretized 2 + 1 dimensional Euclidean space-time. We 
consider a three-dimensional cubic lattice with sites {r} and links {r, //} (/x = 1,2,3 where 
1 and 2 are the space directions and 3 is (imaginary) time). We represent the worldlines of 
the charged particles by means of an integer valued conserved vector field i^f) defined on 
the links of this lattice. Conservation means, naturally, that A M £ /i (r) = at every site r, 
where A^i^ is the lattice divergence. The current loops can be either the worldlines of first 
quantized particles relative to the empty state or the worldlines of excitations relative to a 
reference ground state. Here we envision several physically distinct situations. One system 
we want to describe cosnsits of non-relativistic charged particles with some statistics and 
at finite density. In this case, the worldlines do not end and only close around the time di- 
mension. Another set of cases of interest are relativistic systems (at zero density). In these 
cases, closed loops are allowed and represent particle-antiparticle pairs. Relativistic field 
theories are Charge Conjugation invariant. Non-relativistic systems with or without Charge 
Conjugation symmetry fit this description as well if one now thinks of the loops as represent- 
ing quasiparticle-quasihole pairs. When we compute the electromagnetic response of these 
systems, we will have to bear in mind that, in general, there is a non-trivial contribution 
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from the reference vacuum which must be added to the contributions of the excitations. 

Thus, we are interested in the behavior of theories which have both closed loops and 
with loops that close around the time dimension. The loops can be described in terms of a 
set of bare particles which can be taken to be bosons but their actual statistics is determined 
by a topological term in the action. Typically, this is accomplished by coupling matter to a 
Chern-Simons gauge field 0. Since we want to make all symmetries explicit we choose to 
define the theory on a lattice. Chern-Simons theories have been defined on lattices but they 
are rather cumbersome to work with. Thus, we use an equivalent construction in which we 
directly attach fluxes to particles and, in this way achieve statistical transmutation ]17| . The 
simplest way to incorporate statistical transmutation in a discretized picture is to modify 
the weight of the path integral by a factor which is a pure phase equal to the statistical angle 
S times the linking number of the worldlines of the particles. Theories with the property 
of charged excitations always being bound to magnetic charge (or flux) are said to exhibit 
oblique confinement ||33|| . The Landau- Ginzburg theories of the FQHE |7j are simplest (and 
best understood) example of oblique confinement . 



In order to keep track of particles and fluxes, and to avoid ambiguities, we define the 
worldlines of the fluxes in such a way that they strictly follow the matter worldlines but exist 
on neighboring sites of the dual lattice, which we denote by {R}. Thus, the configuration 
of fluxes on the links of the dual lattice is the same as the configuration of matter on the 
links of the direct lattice, i. e. £^(r) = i^R) where R = f + (|, ~, |). Thus, the vector field 
in represents both the world lines of the particles and the worldlines of the fluxes which, 
therefore, do not have an independent existence. As a consequence, charges cannot detect 
the constituent fluxes of other particles without detecting their charge. 

Notice that, if the fluxes are attached to the worldlines of the particles by a rigid trans- 
lation, then the system is defined in such a way that the particle cannot feel its own flux. In 
this case there will be no phase factor associated with the writhing number (or selfwinding) 
of the worldline. Hence, there is no fractional spin associated with the configurations but 
there will be fractional statistics. However, one can imagine other definitions of the system 
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at short distances. For instance, it is possible to construct the flux-charge composite in such 
a way that the flux rotates around the worldline as the latter winds around itself. In this 
case the flux-charge composite becomes the standard rubber band construction associated 
with the framing of knots. It is well known from the work of Polyakov JZIJ and Witten [19| 



that this construction leads to the conventional spin-statistics connection. With our explicit 
construction there is no fractional spin, so the theory has the invariance 9 — *> —9 discussed 
in the Introduction. We will show below that the presence of this symmetry imposes severe 
constraints on the behavior of the correlation functions. Conversely, if a fractional spin 
is associated with the flux-charge composite, this symmetry will be broken by an amount 
determined by the fractional spin. (We return to this point below.) 
The action S[iJ\ for the field C^if) will be chosen to be of the form 

S[Q = \ E W G ^(r~ W) + \ E W K,Ar, R) UR) 

f,f' r,R 

+ i J2 e(r- i^r) A^f ) + E - #) UR)B^{R') 

R,R> 

+ ^E ■ (2-1) 

The partition function is 

^Enwwi^ ( 2 - 2 ) 

{« f 

where {r*} are the sites of the direct lattice while {R} are the sites of the dual lattice, 
i. e. R — r + (§,§,|)- Here A^{f) is an external electromagnetic field and B^(R) is the 
corresponding field strength. A^{r) couples to matter through two coupling constants, an 
electric charge e and a magnetic charge h. 

The physical interpretation of each term of the action of Eq. (|2.1| ) is as follows: 



1. The matter interactions are represnted by G^ u (f— r 1 ). This is the parity even part 
of the action and it describes the repulsive forces among the particles. For instance, 
the short range part of G flv (f f") essentially suppresses configurations with > 1 
which means that it penalizes the occupancy of one link by more than one particle. 
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Physically, G^if — f") is the propagator of some field that mediates the interactions 
between the particles. The explicit form of the kernel G^if— r*) will be given below. 

2. The statistical interactions induced by Chern-Simons gauge fields are represented by 
the tensor K^if, R!) and it is parity odd. The precise form on a lattice of this tensor 
is rather complicated (see for instance reference fl3"5|). For our purposes, it will be 
sufficient to require that it assigns the correct phase factors depending on the braid- 
ing of the worldlines. The continuum Chern-Simons action (for a gauge field A^) is 
Scsi-Afj) = ^e^uxAfj,d u Ax, where 9 is the Chern-Simons coupling constant. Then |T9|j20[| 
the statistical angle 6 and the coupling constant 9 are related by 5 = 1/(29). The 
weight for a configuration is exp(z<5z^), where is the linking number of the con- 
figuration pC,2l|. Using a magnetostatic analogy, it is possible to find [20] a simple 



formula which relates the linking number ul to a set of currents which describe the 
worldlines. Explicitly one finds [20| 



vl = J d 3 x J d 3 y J^e^x G (x - y) d x J v {y) , (2.3) 

where Gq(x — y) is the three-dimensional green's function, 

- V 2 G (x-y) = 5(x-y) . (2.4) 

Thus we define K^ u so that, to each configuration of matter, the weight of the path- 
integral has a phase factor with argument equal to the linking number of the config- 
uration up to a coupling constant which we identify with the statistical angle 5. In 
other words, the term involving the tensor K^ u represents fractional statistics. By 
construction, it is automatically invariant under periodic shifts of the statistical angle. 
The explicit form of K^ v depends on details of the discretization procedure. We will 
use an explicit form which works provided that the arguments of K^ v are defined to be 
one on the direct lattice and the other on the dual lattice. In continuum description, 
a point-splitting regularization is needed to define the linking number. The lattice 
regularization takes care of this problems automatically. 
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3. The terms in the second line of the action of eq. [2.i| represent the coupling to an external 
vector field A^ir). The space components, 1 and 2, of this vector field can be used to 
implement an external magnetic field. The time component (3) is a scalar potential and 
a constant piece represents a chemical potential (up to an analytic continuation) since 
it counts the net charge flowing forwards in time at a given time slice. In general, A^(r) 
will be an external electromagnetic field that we will use to probe the system. The 
electric coupling e(f) is the conventional lattice version of the usual current coupling 
which we will assume to be mildlyly non-local. The coupling of the vortex current to 
the field strength h(R), is natural as it represents the coupling between the field 
and and the circulating current associated with the vorticity (in the sense of Ampere's 
Law). In any event, this term gets generated automatically under duality even if it was 
not present in the starting hamiltonian. Similarly, we have allowed these couplings to 
be non-local since they will generally become non-local under duality transformations. 
In any case at the fixed points that we will find the non locality of the electric coupling 
is irrelevant. 

4. The last term of the action represents the electromagnetic response of background 
degrees of freedom which do not participate of the dynamics but which affect the 
electrodynamics of the system (e. g., the conductances, etc.) 

Since we are interested in configurations of loops (worldlines) which are large and well 
separted from each other, we do not need to specify the details of the construction of the 
kernels and K^ v at short distances. The precise form of these kernels at distance scales 
of the order of the lattice constant will not affect the universal, long distance, properties. 
Thus, in practice, we will only specify their behavior at long distances or small momenta. 
The only assumption we will make here is that the models of interest can be defined on a 
lattice, in such a way that all the symmetries are preserved. 
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III. DUALITY 



The strategy of this paper is thus to classify the stable phases and phase transitions of 
theories with the action of eq. Q2.1J) (and a number of simple generalizations). The most 



direct way to do this is to use duality transformations [29]. Duality has been used with 
great success in the study of theories with a structure similar to the one we discussed here 
PIj^P%| . Of particular importance to the problems discussed here, is the work of Cardy and 
Rabinovici |H| on theories of two-dimensional statistical mechanics with oblique confine- 
ment. Cardy and Rabinovici showed that, in the context of their models, the phase diagram 
could be derived by means of a sequence of transformations: duality, charge conjugation 
and parity. They were the first to show that, together, these three transformations gener- 
ated the group SL(2, Z) and that the phases are characterized by condensates with oblique 
confinement. The phase diagram of the FQHE exhibits, phenomeno logically, a very similar 
structure. Although a microscopic model to motivate such constructions for systems related 
to the FQHE has not been available, the analysis of the phase diagram by Kivelson, Lee and 
Zhang || and by Ross and Liitken [JTTJ| has followed a line very close to that of Cardy and 
Rabinovici |39| . 



We now proceed to derive a duality transformation of this theory by standard meth- 
ods P5] , pB| -P^]. This amounts to rewriting the theory in terms of a dual theory defined in 
terms of loops on the dual lattice. Since the flux loops reside on the links of the dual lattice, 
duality, as usual, is the same as to exchange electric currents with magnetic fluxes. By 
means of the Poisson Summation formula, the partition function can be written in the form 

{m M } f 

where m^(r) takes values on the integers and £^{r) are now real numbers. 

We will assume that the system obeys periodic boundary conditions in the time direction 
(xs). For simplicity, we will also assume that the system obeys periodic boundary conditions 
in space (this choice does not affect the physics of this problem in any significant way). The 
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configurations of loops can be classified by the number of times they wind around the time 
direction (£3). Since we will always work in the thermodynamic limit and will only be 
interested in systems at zero total center of mass momentum and current, it is not necessary 
to consider configurations which wind around the space directions. However, if we want to 
describe systems at non-zero density (i. e. without Charge Conjugation invariance) we need 
to consider configurations with non-zero winding number around the time direction. 

The constraint on the field t(i{f) can be solved in terms of a topologically trivial piece and 
a field ^ op (r) which accounts for the the loops with non-zero winding number. The topo- 
logical loops are important to describe systems at non-zero matter density . The constraint 
is solved by 

£,(f) = e, UCT A u MR)+^T^ ■ ( 3 - 2 ) 
In what follows, the topological loops £^ op will be set to zero. Notice that, in solving the 
constraint, a gauge invariance has now appeared; since if 4> a {R) solves the constraint, so 
does 4> a (R) + A a F(R), where F(R) is an arbitrary function of R. This gauge symmetry of 
the dual theory will require the use of a gauge fixing condition when we integrated out the 
variables <$> a . 

We now carry out the duality transformation. It is convenient to do it in Fourier space. 
In terms of the Fourier transforms of £^ and M (with wavevector k^), we get 

e ii (k)=ie IH , ff A v (k)<j> < ,(k) , (3.3) 

where 

A v {k) = -i{e ik » - 1) (3.4) 

is a lattice derivative. From now on we will use the notation A u {k) = k u . 

We now specify the kernels G^ v and K^. Without any loss of generality, the Fourier 
transforms of the kernels, G^k) and K^k) can be taken to be of the form 

G»x{k)=g{k)T^{k) 

Kx(k) = M C , x (k) , (3.5) 
15 



where T^xik) and C^xik) are the parity-even and parity-odd ( Cher n- Simons) tensors, 

T^{k) = S^x - C^x(k) = ie n\ P ^7= • (3.6) 



The coupling functions g(k) and f(k) define the theory. Notice that since i^k) is a conserved 
current, a term in proportional to the k^k v does not actually contribute to the action. 

We will consider two models, defined by two distinct choices for g(k). In the first case, 
which we will call Model I, we choose g(k) = g. This choice describes systems with short 
range interactions. In the other case, which we call Model II, we choose G^k) to be long 
ranged, i. e. linn^^o \/k?g(k) = g. As explained in the Introduction, this choice corresponds 
to the non-local (retarded) l/|x| 2 interactions induced by the quantum fluctuations of the 
electromagnetic field (in the full 3 + 1-dimensional space-time) on the particles that move 
on the plane. We will show below that Model II is actually a self-dual theory. The coupling 
constants of these theories will be given by the zero momentum limits of the functions. 

We now complete the duality transformation. The dual theory is defined in terms of the 
integer- valued dual field L^R) 

L^R) = e^ a A u m a (f) . (3.7) 

The dual theory is found by integrating out the field </> M . Because this theory is gauge 
invariant, we must add a gauge fixing term in order to obtain a finite result. The final 
answer is, naturally, independent of the choice of gauge PU . 



Explicitly, we find that the partition function can be written in the form 

Z D [A,(f)] =Z Y: e- s ^ J[6(A,L,(R)) . (3.8) 

{L»} R 

S-ol-Lfj] is the dual action and it is given by 

S dual [L,} = ~J2 L^{R) (G^(R-R') +iK%{R,r>)) L V {R!) 

R,R' 

+ % J2 e D (R - r) L^R) Atf) +J2h D (R- R!) L^R) B^(R') 

+ MVKPrfM?) • (3.9) 
j? pi 
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The dual kernels and charges of eq. (|3.9| ) are in general functions of distance and their Fourier 
transforms are functions of the momentum. We will simplify our notation by dropping, from 
now on, any explicit reference to the momentum dependence of the dual kernels and charges. 
The dual kernels are given by 

Gg = (2n) 2 g D T, v 

K^ = -(2nf^=C^ , (3.10) 
v k 



where the dual coupling functions #d and /d are given by 



9d = (2?r' 



2 



/d = -w'tm^v ' (3 ' n) 



Eq. (|3.11|) is the transformation law of the coupling functions under Duality. Similarly, the 



dual electric cd and magnetic charges hn are 

= — (ef D + k 2 hg 

S.TT x 



„ , -jD 

In 

h D = ^(-eg D + hf D ) . (3.12) 



is the polarization tensor of the dual theory in regimes in which the dual loops 



are suppressed. Thus, it is an effective renormalization of the tensor 11° . In general, the 
full polarization tensor includes contributions coming from both the loops and from the 
background. UP y has the form 

n^ = u^T lv -m^c lv , (3.i3) 

where and 11^ are the amplitudes for the parity-even and parity-odd pieces of H? u , 
11 D -11 + k {e hk j (2 ^ 2 k 2he^ )2 

I#> = - v^(e 2 - h^)^ - Meh^VV (3.14) 

Hence, the polarization tensor acquires additive corrections under a duality transformation. 
This feature will play a very important role in the analysis of the conductance. 
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IV. ELECTROMAGNETIC RESPONSE 



Using the results of the previous section we can, in principle, determine the long-distance 
properties of the possible fixed points. We will parametrize the fixed points in terms of 
their response to electromagnetic perturbations, i. e. the longitudinal conductivity a xx and 
the Hall conductivity a xy for metallic phases and the dielectric constant e and the Hall 
conductivity a xy for the insulating states. 

As usual, these responses are determined from the current-current correlation functions. 
However, the physics of this system is determined naturally by the loop-loop correlation 
functions 

V flu (r-r , ) = (£,(r)£,(7)) . (4.1) 

While is related to the current-current correlation functions, it is not equal to them. 
The current correlation functions are instead equal to the full polarization tensor for the 
system Tl^if — f) which is defined by 



2 



iWf - f) = WO U*)) = - SA J SAJf) iog*«A.(0» 



(4.2) 

A u =0 



where J M is the electromagnetic current. 

Since the current loops {i^} are conserved, i. e. A j[i £ /i (r) = at every site r, we can 
always write T>^ v in terms of transverse tensors. Thus, in Fourier space we write 

2V(*) = V s (k; g, f) T^(k) - iD A (k\ g, f) C^k) . (4.3) 

Here, V s (k;g,f) and T> A {k\ g, f) are the parity-even and parity-odd amplitudes of V^, 
T fa ,(k) and C^ v {k) are the tensors defined in eq. ( |3.6|) . Once again, it should be kept in 
mind that the kernels and the couplings have an explicit momentum dependence. 

The relation between the loop correlation function T>^ v {k) and the current correlation 
function Tl^ik) can be determined by differentiating the partition function. If we also 
separate the U^ u into parity-even and parity-odd pieces, 
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U^(k) = U s (k; g, f) T^(k) - iH A (k; g, /) C^(k) , 



we find 



(4.4) 



U s \k- g, f) = U s (k) - (e 2 - k 2 h 2 ) V s (k; g, f)-2eh Vk 2 V A (k; g, f) 
U\k- g, f) = U A (k) - (e 2 - k 2 h 2 ) V A (k; g, f) + 2eh VtfV s (k; g, f) . 



(4.5) 



A Parity transformation is equivalent to the mapping g — > g and / — > — /. Hence, under 
Parity, the amplitudes V s (k; g, f) and V A {k; g, f) transform like 



V s (k;g,-f)=V s (k;g,f) 
V A (k;g,-f) = -V A (k;gJ) . 



(4.6) 



Thus, the amplitudes II s (k; g, f) and Tl A (k; g, f) transform in exactly the same way. 

By using the definition of the conductivities as linear response coefficients, the compo- 
nents of the conductivity tensor = 1,2) can now be determined by the standard 
relations (in imaginary time x 3 ) 



k 3 ^0 fc ^ K 3 



(4.7) 



Hence, 



a xx = lim — n s (0,/c 3 ) 



'xy 



\hn—U A (0,k 3 ) , 

k 3^° |^3 | 



(4.8) 



where we have set cr n = a 2 2 = <*xx an d (T i2 = _ °2i = &xy We can now relate a xx and a 



'xy 



to the coupling constants of this system 



a xx = a° xx + lim - [-(e 2 - k\h 2 ) V s {k 3 - g, f)-2eh \k 3 \V A (k 3 ; g, f) 
a xy = a° xy + hm \- \-{e 2 - k\h 2 ) V A (k 3 ; g, f) + 2eh \k 3 \V s (k 3 ; g, f) 

K3—>U K 3 L 



(4.9) 



If the system is in a metallic phase both conductivities a xx and a xy are in general different 
from zero. If the system has an energy gap and, hence, it is incompressible, the longitudinal 
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conductivity a xx vanishes. Instead, H s (k)/k 2 has a finite limit as k — > which can be 
identified with the inverse of the effective dielectric constant, e~ l = linifc^ n 5 (/c)/fc 2 of the 
system in an insulating phase. If the insulator has broken Time Reversal symmetry, as in 
the FQHE, the Hall conductivity a xy may be non-zero. 

The phase diagram of the system is spanned by the asymptotic behavior at small mo- 
menta of the coupling functions g(k) and f(k), and of the electric and magnetic charges e(k) 
and h(k). In the next section we will discuss two distinct limiting behaviors for g which 
yield quite different physics. In general, the conductivity tensor at momentum scale k is a 
function of g, f, e and h. In the next section we will exploit the existence of these duality 
transformations to determine the properties of both the stable phases of this system and of 
the phase transitions which separate the phases. 



V. PHASE DIAGRAM 



In this section we will apply the general duality results of section [H| to Models I and II 
and use it to determine the qualitative features of the phase diagram. In this paper we will 
considering only theories with Charge Conjugation invariance. In this case the total charge 
Q is set to zero and only loops that close on themselves are taken into account. 



A. Short Range Interactions 

Let us consider first the case of Model I. In this model, the interactions are short ranged 
and the effective action is characterized by finite limiting values of the coupling constants. 
Hence, we choose the parametrization 

r 9(h) 
q = hm 

fc^o 2tt 

/ = ! im o^? W 

This choice of parametrization reflects the fact that since the interaction is short ranged, 
(7(0)7^0. Also, /(0) is the statistical angle and as such it is defined modulo 2tt. Thus, 
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Periodicity is the shift of / by integers, i. e. f — > / + n. 

Under Duality, the coupling constants of Model I transform like 



9 

9d 



P + k 2 g 2 

/d = ~f 2 + k 2 g 2 ^ 
Eq. ( |5.2[) shows that if f ^ 0, the effective dual couplings as fc — > become 



3d - -~ 



g_ 

f 2 

/d = -j (5.3) 
The electric and magnetic charges obey similar transformation laws 
e D = ^- Jim (ef D (k) + k 2 hg D (k)) = ef D 

/id = ^rlim (-eg D (k) + hf D (k)) = -eg D + hf D (5.4) 

Thus the theory is self-dual in the sense that the original theory and its dual have the 
same form. Clearly, since by shifts of the stastitical angle it is always possible to bring / to 
the fundamental period < / < 1, the effective dual coupling satisfies the inequality gn> g 
and the dual coupling constant is larger that the original one. One usually expects that, if 
the original coupling constant is large, the dual coupling constant should be small; if this 
property holds, the strong and weak coupling regimes are related by duality. In this case, 
the analysis of these limits becomes very simple since in one regime the loops are suppressed 
while in the other the dual loops are suppressed. However, in general, the coupling constants 
of Model I, defined in eq. (|5.2j ) and (|5.3|) , do not have this property. Hence, if the coupling 
constant of the original theory is of order one, the dual coupling is generally larger than 
one. Therfore, for Model I duality does not help (in general) to determine the behavior of 
the system in asymptotic regimes. 

In addition, if / = (up to integer shifts), the situation is quite different. Now the dual 
coupling is long ranged g-o{k) = l/(k 2 g). This behavior leads to logarithmic interactions 
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of matter loops in real space. The fact that a theory of bosonic loops with short ranged 
interactions is dual to a theory of loops with Biot-Savart-like logarithmic interactions is 
well known from the three-dimensional classical XF-model . The loops with long range 



logarithmic interactions are the vortices of the ordered (or superffuid) phase of the XY 
model. 

B. Self-Dual Actions 

The transformation properties of Model II under duality can be best studied in terms of 
the dimensionless coupling constants / (the same as in eq. (|5.1|)) and g which is given by 



f ^ m m ssiim ^m . (5.5) 

For space-time separations x — x', the interaction kernel G^x — x') for Model II Model 
II is transverse and it decays like g/(x — x') 2 . This is the propagator for Maxwell photons 
in the full 3 + 1-dimensional space-time. Thus, Model II describes a system of charged 
particles with fractional statistics moving on a plane interacting with the full quantized 



electromagnetic field (see however, the caveats of reference [28 
For Model II, the dual coupling constants /□ and gr> are 

f 2 + 9 2 f 2 + 9 2 

For the case of Model II, the electric and magnetic charges are required to obey 

e = lim e(k) h = lim \fk^ 2 h{k) . (5.7) 

The dual electric and magnetic charges tr> and ho are 

eD = ef D + hg D 

h D = -eg D + hf D . (5.8) 

It is convenient to define the complex variable z = f + ig. Clearly, while / can take any 
value (although only rational values are physically meanigful), g has to be a positive real 
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number for the interactions to be repulsive and the system to be stable. Hence, the complex 
variable z takes values on H, the upper half of the complex plane. 
In terms of z, Duality is the mapping 

1 . . 

z -> z D = — , (5.9) 

z 

In this notation, Periodicity becomes the mapping 

z — > zb = z + n , (5.10) 

and Charge Conjugation is now 

z^z D = -z* . (5.11) 

Notice that it is always possible to use Periodicity to make Re z = f > 0. Thus, Charge 
Conjugation is actually redundant. 

Hence, Model II is self-dual in the sense that under Duality it maps into a model with 
interactions of the same form. Moreover, it has the extended invariance parametrized by 
the mappings of eq. ( |5.9|) , ( |5.10|) and (|5.11|) . These mappings generate the Modular Group 



SL(2,Z), the group of fractional linear transformations T with integer coefficients of the 

upper half of the complex plane 7i into itself of the form 

_ az + b , 
Tz = — —- 5.12 

cz + d v ' 

with a,b,c,d G Z and ad — be = 1. Clearly, the partition function of Model II is Modular 
Invariant and hence its phase diagram is generated by SL(2, Z). Thus, we expect that 
Model II should have the conjectured structure of the phase diagram of the FQHE PJITf. 
We will show below that Modular Invariance has profound consequences. In particular , it 
determines the possible behaviors of physical observables such as the conductances. 

We also give the transformation rules of the background conductivities under duality 
transformations. They are 

1 



D _ 

^ xx ' 

D _ , J_ 

°xy °xy 27]- 



[e z - h z )g D - 2ehf D 
-(e 2 -h 2 )f D -2ehg D \ . (5.13) 
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It is convenient to rescale the amplitudes of Eq. ( fl9|) and to define 

?ir 9-7T 

d s (g,f) = ^=V s (g,f) d A (g,f) = ^=V A (g,f) . (5.14) 



The self duality of the theory imposes constraints on rescaled amplitudes d s (g,f) and 
d A (g,f) and, through them, on the physical conductivities. Thus, from Eq. (f4.9|) , we get 
that the physical conductivities a xx and a xy at coupling constants g and /, for specified 
electric and magnetic charges e and h are 

1 



-(e 2 -h')d b (gJ)-2ehd A (gJ) 
_( e 2 -h 2 )d A (gJ) + 2ehd s (g,f)] . (5.15) 



Likewise, for a Time Reversal Invariant system, we have the further restriction imposed by 
Eq. (4J3) on the amplitudes V s and V A , 



d s (g,-f) = d s (gj) 

d A (g,-f) = -d A (gJ) . (5.16) 

We stress here that Time Reversal follows from choosing a specific definition of the loops 
at short distances in such a way that no fractional spin is assigned to a system with a 
given fractional statistics. For systems with fractional spin, Time Reversal is broken and the 
Reflection Symmetry of Eq. ( |5.1(j| ) does not hold. 

The transformation properties of the quantities of physical interest become very transpar- 
ent by using complex notation. Thus we define the complex conductivity S, the amplitude 
D and the charge H, 

£0) = o- xy (f, g) + ia xx (f, g) 
D(z) = d A (f,g) - id s (f,g) 

H = h-ie , (5.17) 

where, for simplicity, we have used the notation D(z) = D(z,z*). ( Here the stars denote 
complex conjugation.) Naturally, this notation does not imply that D(z) is necessarily only 
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a function of z and thus that it is analytic. In general, the amplitude D is a separate function 
of z and z* since, a priori, there is no reason to expect D to be an analytic function of the 
coupling constants / and g. However, as we will see below, there are situations of physical 
interest where it can be locally analytic. We will also see that analyticity has far reaching 
consequences. In what follows we will write all functions of / and g as functions of z. (The 
same caveats on analyticity still apply.) 

With this notation, the complex conductivity £ becomes 

£(z)=E + ^ir(z) . (5.18) 

Similarly, we find that the duality transformation for the charges is now simply 

H D = Hz* D . (5.19) 

Since this is a self-dual thery, the conductivity of the theory at z and at Zd must be the 
same. Hence 

= S D + ^D*(z D ) , (5.20) 

where 

H 2 

Z D = Z 0+ Z * D . (5.21) 

The consequence of these relations is the the amplitude D(z,z*) obeys the following func- 
tional equation 

z 2 D{z) + z = £>(--) , (5.22) 

z 

which is the transformation law of the amplitude D(z) under duality. In addition, D(z) 
must periodic which here means 

D(z) = D(z + 1) . (5.23) 

Similarly, the Charge Conjugation symmetry transformation of the amplitude under 
{fi9) ~ * {~fi9)i Eq. ( |5.16| ), which only holds if the system is invariant under Time Re- 
versal, becomes now the Reflection symmetry, 
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D(-z*) = -D*(z) , (5.24) 

which expresses the invariance under Time Reversal. Clearly, D(z) cannot be analytic if the 
system is Time Reversal Invariant. 

In the next subsection we will use Modular Invariance (and Time Reversal when appro- 
priate) to determine constraints on the allowed values of D(z). Once D(z) is known, the 
conductance £(2) is determined uniquely by Eq. ( |5.21| ) once a choice of the electromagnetic 
coupling H = h — ie is made. In what follows, we will choose the bare model to have H = —ie 
(i. e. h = 0). 



C. Consequences of Modular Invariance 

Let us now discuss the consequences of Modular Invariance. The Modular Group 
SL(2, Z) is an infinite, discrete, non-abelian group. It has two generators, which we will 
take to be the Periodicity P and the Duality S. Please notice that in the mathematics 



literature [}42j] , as well as in String Theory the convention is to call T-duality what we 
call Periodicity and 5-duality what we call Duality. By the action of SL(2, Z), it is always 
possible to map an arbitrary point of the upper half-plane 7i into the Fundamental Domain 
F = {z E H\\z\ > 1, |Rez| < |}. Particular cases include the mapping the interior of the 
unit circle \z\ < 1 to the exterior \z\ > 1, which is the conventional mapping of weak and 
strong coupling regimes into each other. 

The unit circle C\ = {z; \z\ = 1} plays a special role since it is a self-dual manifold 
although it is not modular invariant. On the unit circle f 2 + g 2 = 1, duality reduces to 
the mapping fp = —f and = g. In particular, the self-dual boson point on the unit 
circle is (/, g) = (0,1) = i and the self-dual fermion point on the unit circle is (/, g) = 
(1/2, y/3/2) = p. The boson and fermion self-dual points have very special properties. Let 
[i] and [p] denote the sets of points of the upper half-plane TC which are the images of % and p 
under SL(2, Z) respectively. For similar reasons, we will call the sets of points with integer 
real parts ( i. e. f £ Z) the boson lines. Likewise, the the sets of points with half-integer 
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real parts i. e. f E Z + |) are the fermion lines. In particular, Ke[z\ = 1 will be the boson 
line while Re[z] = |is the fermion line. 

In addition of the finite modular fixed points [i] and [p] the modular group also has fixed 
points at extreme regimes of the upper half plane. The extreme fixed points are the point 
{ioo} and [oo],the set Q of rational points on the real axis. The point {ioo} is the limit 
Imz — > oo and corresponds to the strong coupling limit g — > oo. In this regime the loops 
are suppressed. Conversely, [oo] is the set of points of H with g — and / rational. These 
points map into each other under SL(2, Z) and are also images of the point at infinity on 
the real axis. Also, since g = 0, physically the set of points [oo] corresponds to systems of 
non-interacting charged particles. Also, since g — > 0, for these points it is not legitimate to 
neglect other, shorter range, interactions. In fact, close enough to these points, Model II 
behaves like Model I. Clearly, the behavior near these points is quite simple since, for large 
values of g the loops are suppressed while for small values of g the dual loops are suppressed. 

The bosonic point i and the fermionic point p are examples of fixed points of SL(2, Z) 
and describe interacting systems. It can be shown that the set C of all points of the 



upper half-plane TC which are fixed points under the Modular Group SL(2, Z) (namely the 
set of points in 7i which are fixed points of a modular transformation), is the union of [i], 
[p] and [oo], i. e. C = [i] U [p] U [oo]. In what follows we will say that a fixed point is bosonic 
or fermionic if it is in [i] or [p] respectively. Notice that there will be points on fermion lines 
which are bosonic since they are images of i, for example z — | + |. 

These three families of modular fixed points should determine the qualitative global 
features of the phase diagram. This fact has been stressed by Ross and Lutken [jTlj] , and 



by Lutken |j| who have given a detailed description of these modular fixed points. In 



that work, they proposed a set of phenomenological Renormalization Group flows and fixed 
points which are compatible with the requirement of modular invariance. Thus, Ross and 
Lutken were led to the conclusion that the RG fixed points should also be SL(2, Z) or 
modular fixed points. Thus, the set [oo] describes plateaus while the finite fixed points [i] 
and [p] describe phase transitions. The qualitative phase diagram of Model II has the same 
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qualitative structure as the phase diagram of Liitken ]12| (see in particular Liitken's figure 
2). This is also the phase diagram that has been proposed by KLZ for the Fractional Quantm 
Hall Effect ||. However, we will show evidence that suggests that Model II has properties 
that we would normally assign to critical points at most points of the parameter space TC. 

Before proceeding further it is convenient write down the transformation law of the 
complex conductance D(z) for an arbitary modular transformation. The transformation 
laws for D(z) under duality Eq. (|5.22|) and periodicity Eq. (|5.23|) , generalize for an arbitrary 
modular transformation T G SL(2, Z) of the form of Eq. ( |5.12| ), to the following law 

D(^-^) = (cz + d) 2 D(z) + c(cz + d) . (5.25) 

Once again, we stress that so far we have not made any assumptions on the analytic prop- 
erties of D(z). 



Equation ( |5.25|) has a very striking and suggestive form. Since the partition function 
is invariant under the action of SL(2, Z), naively we would have expected that the complex 
amplitude D(z) should also be invariant under SL(2,Z). However, the invariance of the 
partition function does not require the invariance of correlation functions. Instead, only 
covariance is actually required. This is precisely what happens in critical phenomana where 
scale invariance only requires that correlation functions to exhibit scaling behavior. Thus, in 
principle, we could have anticipated that if D(z) was not a modular invariant, it would in- 
stead transform homogeneously under the action of the Modular Group SL(2, Z). Functions 
F k k (z, z*), labelled by two integers k and k, that transform homogeneously under SL(2, Z) 



as 



F k> - k (Tz, (Tz)*) = (cz + d) k (c*z* + d*) k F krk (z, z*) (5.26) 

are said to be modular of weight (k,k). In the special case when the function F(z) is 
holomorphic (and k — k), it transforms like 

F 2k (Tz) = (cz + d) 2k F 2k (z) (5.27) 
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and it is called a modular form of weight 2k. Typical examples of modular forms are the 
Eisenstein series G 2 k( z ) 

G ^ z ) = E !_ , 2k ( 5 - 28 ) 

[nz + mY K 



and their generalizations f44 



However, Eq. (|5.25|) transforms inhomogeneously, i. e. in addition to the scale factor of 
the form of Eq. ( fT27D (with weight 2k = 2), there is an extra additive contribution. This 
feature can be traced back to the duality transformation properties of the polarization tensor. 
Indeed, as Eq. (|3.14j ) shows, an external electromagnetic perturbation (no matter how weak) 
breaks self-duality. We will see shortly that D(z) is related to a (loosely defined) modular 
form of weight 2. This inhomogeneous transformation law has a striking similarity with the 
transformation law of the energy-momentum (or stress-energy) tensor in two-dimensional 
Conformal Field Theory which exhibits a Conformal Anomaly |¥5 |. By analogy we can 



regard the inhomogeneous transformation law of Eq. ( |5.25| ) as a Modular Anomaly. 

It is straightforward to show that if zq is a fixed point, the generalized transformation 
law Eq. ( |5.25| ) determines uniquely the value of D(z ) at the fixed points. Indeed, let z e C 
and let T zo be the modular transformation that leaves z invariant, i. e. T zo z = zq. Then 
Eq. ( |5.25|) reads 



D(T Zo z ) = D{z ) = {cz + d) 2 D(z ) + c(cz + d) . (5.29) 

Hence we find 

= l-iczo + W ■ (5 ' 30) 
The value of D(z ) of Eq. ( |5.30| ) determines uniquely the allowed values of the conductivities 
at the fixed points. However, it is straightforward to check that 

D(z) = — — (5.31) 
v ; 2lmz v ; 

satisfies both the requirement of periodicity P and it transforms correctly under duality S. 

Hence, the value of D(z) at the fixed points given by Eq. ( |5.30|) must be equal to 
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Therfore, D(z) is pure imaginary at the fixed points. 

Let us derive its consequences for the specially important cases of bosons and fermions: 

a ) Bosons: The boson fixed point z = i is left invariant by the Duality transformation, 
which is the SL(2,Z) transformation with a = d = and b = — c = 1. Eq. (|5.30|) 
yields the value 

D(i) = l - (5.33) 

which is pure imaginary. The conductivities for a theory of self-dual bosons become 

e 2 e 2 
o- xx = a° xx - —d s = a° xx + — 

e 2 

°xy = oly - —d A = o-° xy . (5.34) 

Notice that the odd-amplitude d A = 0. Except for the background conductivities a xx 
and a xy , the conductivity a xx of self-dual bosons is indeed universal and equal to 
This is the conjectured value of the universal conductivity. 

Likewise, a generic bosonic fixed point (i. e. a fixed point zq e [i]) has the form 

, _ at + b . , 

z Q = Ti-+ Zo % = . . (5.35) 
ci + a 



Eq. ( |5.32| ) tells us that D(z ) is pure imaginary. The longitudinal conductance at a 
generic boson fixed point is universal although it is not equal to 

Fermions: The fermion fixed point z = p = \+i^ is left invariant by Duality followed 
by Periodicty (with period 1). This is the modular transformation with a = — b = c = 1 
and d = 0. For this point Eq. ( |5.30| ) yields the value 

D(p) = i^- (5.36) 



30 



which, once again, is pure imaginary. The conductivities for self-dual fermions are 
then given by 



^ j ^ 

0~ X X = C M — — US — a rrrr. + ~ 



" 2Tr~° ' 2;-; 3 

e 2 



Vxy = °xy ~ ^ d A = 0%, ■ (5.37) 

Thus, the conductivities are again universal since, up to effects of the background 
conductivites, we get a xy = and a xx = f~ 3 as the value of the critical conductivity 
for this self-dual fermion system. This analysis can be extended to all the fermionic 
fixed points, zo G [p]. Once again Eq. (|5.32|) yields a universal conductace. 

It is interesting to consider the fixed points (both bosonic and fermionic) on the fermion 
line z — I + ig. These points are mapped into other points on the same fermion 
line z' = \ + ^ by the modular transformation (a,b,c,d) = (—1,1,-2,1). This 
transformation has one fixed point on this fermion line, namely the bosonic point 
z = I + The transformation law of the amplitude D(z) now becomes 

D (\ + 7- ) = ~ A 9 2 D C- + ig] + . (5.38) 
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It is easy to check that D f| + |J = i. Likewise, z = p maps onto z = | + i^-. The 
transformation law now yields D { \ + i^) = i^>- Once again, we find a universal 



conductance. 

c ) Extreme Fixed Points: These are the fixed points in the class [00] U {ioo}. It is 
straightforward to show that D(z) vanishes near {ioo}. For instance, for z — f + ig 
and g — > 00, the contribution of the loops to the partition function is suppressed. In 
the absence of loops the loop correlation functions vanish. More specifically, since the 
model is periodic in the statistical angle /, the partition function can be written as a 
sum over linking numbers or, what is the same, as the Fourier series 

Z(J,g)= £ Z UL (g)e 2mfUL (5-39) 
v L ez 
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Here, the integer ul is the linking number and Z„ L (g) is the partition function at fixed 
linking number. Z UL (g) is positive and it is a function of g alone. For large but finite 
g, each term of this sum is of the order of exp(— gSo). The the action Sq depends on 
the shape of the loops but it is finite for large loops. Hence, for large g, only small 
loops will contribute to Z UL (g). In this regime, this series is convergent and, in fact 
the partition function at statistics / is bounded from above by the partition function 
for bosons (which also has an exponential dependence in g, at large g). Intuitively, 
this argument seemingly also implies that the leading non-zero contributions to the 
correlation functions are also exponentially small in g. If this were so, in this regime 
D(z) would have an essential singularity in g, for large g. While this is perfectly 
correct for a system of bosons, it is also conceivable that for generic statistics / the 
phase fluctuations could conspire to modify this behavior even at large g. In fact this 
is what should happen if Eq. Q5.32D would hold even at large values of g and not just 



at the fixed points. Without further elaboration we conclude that, in any case, we 
should expect D(z) to vanish at large Imz, either eponentially fast or with a weaker 
power law dependence. In all scenarios, the physics of this system near {ioo} is clearly 
that of an insulator. 

On the other hand, near the other extreme fixed points [oo] (on the real axis g — ► 0), the 
loops proliferate. However, since near these points the dual loops are suppressed, the 
behavior of the system in these regimes can be determined by using the transformation 
law. For instance, Duality maps the boson points z = ig onto the boson z = —i/g. 
For these points, the transformation law reads 

D{-) = -g 2 D(ig)+ig . (5.40) 
9 

Hence, if g — > 0, this requires the divergent behavior D(ig) « i/g. Similarly, the points 
on the fermion line z = \ + ig for small g are mapped by the modular transformation 
(—1, 1, —2, 1) onto points on the same line but now with large values of g. Finally, 
since in the extreme limit g — > the long range interactions effectively disappear, 
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the behavior of the system is dominated by short range interactions and/or mass 
terms. In particular, the free Dirac fermion point g — > has a Parity Anomaly whose 
precise value is determined by the sign of the (infinitesimal) mass m term . In our 
language, the Parity Anomaly is the finite real limit D(z = \ + i0 + ) = sign(m)|. We 
will not elaborate on this interesting issue here any further. In any event, even along 
the fermion line we expect essential singularities for large g (up to effects of the Parity 
Anomaly). Notice, however that since the modular fixed points accumulate on the real 
axis, there will always be a lenght scale in which the long range interaction dominate 
no matter how small g may be. 

The analysis we have just presented shows that, at least for the fixed points, D(z) is 
completely determined and its values at the fixed points are universal and independent of the 
analytic properties that we may assign to it. Also, and in the same vein as in the discussion 
of the Introduction, the fact that a xx is neither zero nor infinite in this system (which is 
not disordered) implies that, at the modular fixed points, the system must necessarily be 
critical. This argument suggests very strongly that the system must also be scale invariant 
at these fixed points. 

It is natural to ask if these are the only points of the upper-half plane 7i where a xx 7^ 0, 
or, equivalently, ImD(z) 7^ 0. However, now we have to confront the problem that the 



transformation law of Eq. |5.25| alone does not determine completely D(z) any longer. To 
make further progress we need to know something about the analytic properties of D(z). 

We have already argued above that if a fractional spin is not assigned to the particles, 
this system has a time reversal invariance which requires the reflection symmetry of the 
amplitude D(z). In turn, reflection symmetry implies that D(z) cannot be an analytic 
function of z in the sense that it is a function of both z and its complex conjugate z* . It is 
natural to assume that there couldbe a spin-statistics connection for which D(z) is analytic. 
It is also tempting to speculate that D(z) would become analytic precisely for the natural 
spin-statistics connection (although we have no proof of this) . 
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Let us assume momentarily that D(z) is analytic ( holomorphic) and derive the con- 
squences that follow from this assumption. It turns out that, on the upper-half plane 7i, it 
is impossible to construct a holomorphic function of D(z) which transform homogeneously 
under SL(2, Z) like a modular form of weight 2 since this would violate the Residue Theo- 
rem. However, in the theory of modular forms [ P2"| , [P| it is proven that instead there exists 
a unique holomorphic function which, remarkably, transforms precisely as required by our 
inhomogeneous transformation law of Eq. ( |5.25|) . Thus, even though Model II may not 
require analiticity it is certainly compatible with it. The precise statement of this result is 
that D(z) is a holomorphic function of z (hence, only a function of z) of the form 

D(z) = -ATtiG 2 {z) . (5.41) 

Here G 2 {z) is the holomorphic function defined by the Fourier series 

1 oo 

G 2 (z) = -- + Y,vi(n)e 2mnz , (5.42) 

Z4 n=l 

where o~i(n) is the sum of all the divisors of the positive integer n, 

^iW = E r • ( 5 - 43 ) 

r\n 

In terms of the variable q = exp(2niz), the first terms of this Fourier expansion are 

G 2 {z) = -^ + q + 3q 2 + 4q 3 + ... . (5.44) 

Alternatively, G 2 (z) can be written eularized Eiesenstein series of the form 

G 2 {z) = l — lim Y '- — + — — (5.45) 

87r" ! e^o+^ [mz + n) 2 \mz + n\ e 87rlm2; 



The significance of Eq. ( |5.41| ) is that, if analyticity were to hold exactly, the complex con- 



ductance would be a universal function of z — f + ig, completely determined on the entire 
upper-half plane. 

In contrast, if analyticity does not hold, many options are available. Indeed, we have 
already noticed that 
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is a solution. It also happens to satisfy reflection symmetry and it is also consistent with all 
the values of D(z) that we have determined so far, including the asymptotic behavior for 
g — > 0. However, we now have infinitely many solutions available to us. Indeed, let F(z) 
be a modular invariant (not necessarily holomorphic) function of z. An example of such 
functions is the Jacobi Invariant J(z) (which has an essential singularity at infinity). The 
derivative F\z) is modular of weight 2 (but not holomorphic). Then we have 

C « = 2ib + F 'M- < 5 ' 47 > 

For reflection symmetry to hold it is sufficient to write 

fl W = 2lk + !|F ' WI ' (5 - 48) 

which satisfies all the requirements. 

Thus, analyticity is a very strong (and highly desirable!) constraint which determines 
uniquely the amplitude D(z) and the conductances for all points of the upper half-plane. 
Analyticity appears to be related to the existence of a spin-statistics connection in a manner 
which is not yet fully understood. If analyticity holds, our results show that it is possible 
to have a more general (and interesting) type of superuniversality: the conductances are a 
universal function of the coupling constant g and of the statistical angle /. In a sense, the 
complex conductance behaves very much like a universal amplitude in critical phenomena, 
which are universal functions of effective coupling constants. The latter, of course, are related 
in a highly non-universal manner to the coupling constants that define the physical systems 
at microscopic scales. The same structure appears to hold here. We close by stressing that, 
even if analyticity does not hold, modular invariance alone is sufficient to prove that the 
complex conductance is different, but determined, at the different modular fixed points. 
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VI. PERIODICITY IN QFT AND CRITICAL BEHAVIOR 



In previous sections we derived a number of duality transformations for two classes of 
theories of interacting loops endowed with statistical phases. 

We showed that, for Model I and for a subset of possible values of the statistical angle 
/, the theory with short range interactions among loops is dual to a theory with long range, 
logarithmic, interactions among loops. The theory with short range interactions among loops 
is essentially identical to the high temperature phase of the 3D XY classical Model. In the 
high temperature phase, the loops are the loops of the diagrams of the high temperature 
expansion. These diagrams are also a lattice version of the Feynaman diagrams of a theory 
of a self-interacting complex field. Close to the critical coupling constant (or T c ) this theory 
becomes equivalent to |0| 4 theory. This picture is known to hold from detailed studies of 



critical behavior both within the 4 — e expansion and directly in three dimensions [55 



The dual theory is the low temperature phase of the 3D XF-Model. In this regime the 
loops are vortex defects of the order parameter field. The fact that a theory with long 
range interactions among loops is dual to a theory with short range loop interactions, both 
separated by the same second order phase transition, indicates that the long range field of 
the vortices of the low temperature phase is actually screened and that the actual fixed point 
is that of the |0| 4 theory. The critical behavior of this theory is controlled by the Wilson- 
Fisher fixed point of its 4 — e expansion. It is well established that the phase transition of the 
3D XY model is second order. The disordered phase can be viewed as a vortex condensate 
but, unlike the 2D XY model, it is not a vortex-liberation phase transition. In fact, as the 
transition is approached, core size of the core of the vortices grows large and the long range 
interaction of the vortices becomes screened. The transition of the 3D XY model is actually 
controlled by the Wilson-Fisher fixed point of |0| 4 Landau theory. For other values of the 
statistical angle /, the theory is equivalent to a |0| 4 theory coupled to a Chern-Simons gauge 
field with a suitable value of the Chern-Simons coupling constant. 

If matter is second-quantized, the situation is somewhat less clear. For second-quantized 
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theories defined on a lattice [4^- ^9 , 35] , Periodicity always holds and so does Charge Con- 



jugation in the appropriate cases. For non-relativistic dense systems in two-dimensional 
continuum space (i. e. all Chern-Simons theories of the FQHE) are invariant under P and 
invariance under C only happens under special circumstances. For continuum relativistic 
field theories the situation is less clear since the symmetries depend on subtle properties of 
this theory at short distances and are thus strongly sensitive to regularization prescriptions. 
For instance, field theories with a single relativistic fermi field with a four-fermi interac- 
tion and coupled to a Chern-Simons gauge theory have been considered recently by several 
authors [^0|-|52|. Within a perturbative approach, which uses dimensional regularization, 



these authors study the critical behavior of such theories. In particular, it is found |52| that 
the dimension of the four-fermi operator A 4 is a continuous function of the statistical angle 
5 (related to the Chern-Simons coupling constant 9 = 1/25) and that, in particular, A 4 
decreases like A4 — 4 ~ 5 2 (here 4 is the dimension of the four-fermi operator at the free 
field fixed point). Moreover, due to the topological charcter of the Chern-Simons action, 



its coupling constant acquires at most finite renormalizations p3[ and, consequently, it does 



not flow under the renormalization group. Thus, if taken at face value, these results seem 
to imply that the critical properties of at least relativistic matter fields coupled to Chern- 
Simons gauge fields do depend parametrically on the value of the Chern-Simons coupling 
constant 9 and that these transitions are neither necessarily superuniversal nor periodic in 
5. 

The results of references f)0H52| hold in a perturbative expansion in powers of 5 (i. e. the 
deviation from fermions). Hence, periodicity cannot possibly be respected even if the full 
theory had this symmetry. However, the use of dimensional regularization imposes more 
serious limitations. As a regularization procedure, dimensional regularization yields unam- 
biguous results only for theories (and observables) which can be continued in dimension in a 
natural way. However, dimensional regularization is notoriously unreliable for the study of 
theories with operators which involve the Levi-Civita epsilon tensor, such as Chern-Simons 
gauge theory itself and many others. The dimensional continuation of operators with the 
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epsilon tensor is plagued with ambiguities with additional (or evanescent) operators which 
make important contributions even though they formally do not exist in three dimensions. 
From a physical point of view, the vanishing of all divergent contributions other than log- 
arithmic, which in standard problems is an advantage of dimensional regularization, turns 
into a problem in this case. For instance, in these theories there are ultraviolet linearly 
divergent contributions to the fermion self-energy which induce a non-zero mass for the 
fermion. Since fermion masses break parity and time reversal in 2 + 1-dimensions, they alter 
the symmetries of the system. Thus, different choices of cutoff may affect the physics. For 



instance, in reference |54j] the anisotropic quantum Heisenberg antiferromagnet was studied 
using a lattice Chern-Simons approach. This theory is a lattice version of the relativistic 



theory. It was found |5J] that the parity breaking induced fermion mass terms are essen- 
tial for the spectrum of the theory to agree with the spectrum rigurously known for the 
anisotropic quantum Heisenberg antiferromagnet. 

In any event, a full non-perurbative construction of these continuum field theories in 
which Periodicity is an explicit symmetry is still missing. Our discussion in the rest of the 
paper suggests a possible way to do this. Our results indicate that perturbation theory may 
be misleading. 



VII. CONCLUSIONS 

In this paper we constructed of a set of models which qualitatively describe the transi- 
tions between quantum Hall plateaus and related problems. We gave an explicit construction 
of a model (Model II) which is invariant under periodic shifts of the statistical angle and 
under duality transformations. In fact the model is invariant under the full Modular Group 
SL(2, Z). We also investigated the transformation laws under SL(2, Z) of physical correla- 
tion functions. In particular we showed that the complex conductance transforms (loosely 
speaking) like a modular form of weight 2. This observation alone enabled us to determine 
the values of the complex conductance at the fixed points of the modular group. The values 
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of the conductances at these fixed points are universal but not constant. We present an 
argument that shows that the system is critical at least at these points. In a sense the 
conductance behaves like a universal amplitude of a critical system. Further assumptions 
about symmetries have to be made in order to determine the conductance away from the 
fixed points. If the complex conductance is a holomorphic function of z — f + ig, then the 
modular transformation laws determine the conductance completely on the entire upper half 
of the complex plane. This is perhaps the strongest form of superuniversality. In fact, all 
the non-universal physics is encapsulated in the dependence of the (long distance) effective 
coupling constant g on the microscopic coupling constants. For Time Reversal Invariant 
systems we find that symmetry arguments alone are no longer sufficient to determine the 
conductance. It is quite possible that the system should have a hidden supersymmetry for 
D(z) to be holomorphic. 

A number of additional issues remain to be resolved. The models were constructed in 
such a way that particle-hole symmetry is exact. Also, we have assumed that there was 
no external magnetic field (on average). Clearly, in order to describe the full set of phase 
transitions of the FQHE these assumptions need to be removed. In the context of Model I 
this has already been done in all of the studies of the FQHE plateaus. However, for Model 
II we have chosen not to consider this generalization here. Although in principle this is 
straightforward, this generalization already poses an important question. In the context of 
Model II we have shown that at the modular fixed points that complex conductance is pure 
imaginary and this implies the vanishing of the Hall conductance in a modular invariant 
system. This result followed from the fact that the complex conductance transforms like a 
modular form of weight 2. We have also observed that, in the presence of an electromagnetic 
field, the partition function is no longer a modular invariant. At finite fields and chemical 
potential, the partition function will transform in a non-trivial way. Thus, away from perfect 
particle-hole symmetry, the system can have a non zero value of a xy at the fixed points. 
A related issue is the description in this language of the Parity Anomaly in theories of 
relativistic fermions. 
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Finally, there is the question of the unitarity of Model II. Since Model I describes a theory 
of particles with fractional statistics with short range interactions it can also be written a 
local field theory of a complex scalar field coupled to a Chern-Simons gauge field. Although 
the phase transitions of these theories are not understood, their behavior at their stable, 
semiclassical regimes is well known. In particular, for rational values of the statistical 
angle, these theories are unitary. In contrast Model II, which turned out to be a more 
interesting and richer theory, cannot be mapped onto a local field theory for generic values 
of the statistical angle. Nevertheless, on any bosonic line, the partition function has positive 
weights and it is unitary. However, on any fermionic line / = \ + n (with n G Z) Model 
II is also unitary since, on such lines the theory can be related to a non-local version of a 
Maxwell-Chern-Simons theory. Since Model II is modular invariant, it follows that unitarity 
should also hold on the images of all of the bosonic and fermionic lines under the full action 
of SL(2, Z). Therefore, we expect unitarity to hold on a dense set of points of the upper-half 
plane. However, it is quite possible (and likely) that unitarity may not hold for irrational 
values of the statistics. 
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